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Abstract 



A numerical solution of a generalized Swift- Hohenberg equation in two dimensions 
reveals the existence of a spatio-temporal chaotic state comprised of a large number of 
rotating spirals. This state is observed for a reduced Rayleigh number e = 0.25. The 
power spectrum of the state is isotropic, and the spatial correlation function decays 
exponentially, with an estimated decay length ^ 2. 5 Ac, where Ac is the critical 
wavelength near the onset of convection. Our study suggests that this spiral defect 
state occurs for low Prandtl numbers and large aspect ratios. 
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The spatio-temporal chaotic behavior of spatially extended, dissipative systems 
has been intensively studied in recent years The transition to spatio-temporal 
chaos has been observed in various physical systems such as Rayleigh-Benard con- 
vection 0, optical instabilities 0, flames 0], and chemical reactions |Q. Spatio- 
temporal chaos manifests itself through a breakdown of global spatial coherence. 
However, a macroscopic coherence length -a length scale below which the pattern 
appears coherent- may still be observed. This form of chaos is often referred to as 
weak turbulence. Recently, a chaotic spiral defect pattern (which we will refer to 
loosely as spiral chaos below) has been observed in Rayleigh-Benard convection in a 
non-Boussinesq fluid (CO2 gas) 0, for a large aspect ratio and at moderate Rayleigh 
numbers. In previous experiments on convection in CO2 gas |]^, the spontaneous 
formation of a stable rotating spiral pattern was observed at a lower Rayleigh num- 
ber. The chaotic state observed more recently is comprised of a large number of 
rotating spirals. Spirals nucleate, interact and annihilate yielding a macroscopically 
disordered pattern. In this paper, we show that a generalized Swift-Hohenberg equa- 
tion which includes a quadratic nonlinearity and coupling to mean flow can account 
for the formation of this chaotic spiral pattern. We find that chaotic spiral patterns 
are spontaneously formed during the transition from the conduction state to rolls, in 
agreement with the experimental observations. 

It is important to distinguish from the outset among the chaotic dynamics of dissi- 
pative systems that involve a small number of degrees of freedom, the chaotic spatio- 
temporal behavior of systems with many degrees of freedom, and fully developed 
turbulence. We follow Hohenberg and Shraiman who defined three characteristic 
length scales: the dissipation length In, the excitation length Ie and the correlation 
length ^. The ratios of these lengths with respect to each other and with respect to 
the characteristic system length L are used to characterize the chaotic state. The 
dissipation length //) is the characteristic length at which energy is dissipated. The 
excitation length Ie is the characteristic length in which energy is injected into the 
system. Thus, for example, near the onset of Rayleigh-Benard convection in an in- 
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finite system, In ^ k^-^, where kc is the critical wavenumber. The excitation length 
d k'^^, where d is the thickness of the fluid layer. For a small aspect ratio and 
for moderate values of the reduced Rayleigh number R/Rc (where Rc is the critical 
Rayleigh number), we have L^I^^Ieso that only a few modes are excited. Hence 
the dynamics of such systems can then be explained by considering only a small num- 
ber of modes. At large reduced Rayleigh numbers R/Rc » 1, one has Id « ^e, 
which is the fully developed turbulent regime. We focus here on the case of moderate 
R/Rc, but a large system L >> Ie- If ^ > L, we would be still effectively dealing with 
a small system, coherent in space but that may be either regular or chaotic in time. 
When ^ << L, on the other hand, the dynamical behavior is incoherent in space, 
and the system may be chaotic both in time and space. This regime that occurs for 
L » Ie ^ Id ^ d a.t moderate R/Rc, is the regime of spatio-temporal chaos or weak 
turbulence, and is characterized by a dominant macroscopic coherence length ^. 



We model convection in a non-Boussinesq fluid by a two-dimensional generalized 
Swift-Hohenberg model [10|, defined by Eqs. ([l|)-@ below, which we solve by 
numerical integration. The Swift-Hohenberg equation and various generalizations of 
it have proven to be quite successful in explaining many of the features of convective 
flow in fluids, particularly near onset |Tl], |12|, [ 



1^ . As we show in this paper, the 



same holds true for the formation of spiral chaos in non-Boussinesq fluids. Our model 
is defined in dimensionless units by 

d'ip{f, t) 



dt 



V" + l 



where U is the mean flow velocity. 



(3) 



The boundary conditions are. 



2p\B = n- V^ls = C\B = n- VCIb = 0, 
4 



(4) 



where n is the unit normal to the boundary of the domain of integration, B. Equation 
(HI) with g2 = Qm = ^ reduces to the Swift-Hohenberg (SH) equation. The scalar order 
parameter ?/'(r, t) is related to the fluid temperature in the mid-plane of the convective 
cell, and C{r,t) is the vertical vorticity potential. Mean flow arises when the vertical 
vorticity is driven by roll curvature and amplitude modulations. Coupling to mean 
flow has been shown to play a key role, for example, in the onset of weak turbulence 
in Boussinesq fluids [0, |15|, |16[. The quantity e' is the scaled control parameter, 
e' = (fc^)e, where e = — 1) is the reduced Rayleigh number. Here R is the 
Rayleigh number, Rc is the critical Rayleigh number for an infinite system, kc is the 
critical wave number, is a characteristic length scale, Pr is the Prandtl number, 
and is an unknown constant. 



The values of the parameters that enter the equation have be chosen in the range 
appropriate for earlier experiments of Bodenschatz et al on CO2 0- In order to 
estimate them in terms of experimentally measurable quantities, we have derived a 



three mode amplitude equation from the generalized Swift-Hohenberg equation [14 



From the experiments described in |^, we have estimated (72 ~ 0.35 and Qm ~ 50. The 
value of e' used in the numerical solution is 0.7, which is related to the experimental 
value e in Ref. by e = 0.3594e'=0.2516. We have chosen = 2 to simulate 
approximately the experimental rigid-rigid boundary condition. (Note that = 
corresponds to a free- free boundary condition.) We note that for these values of the 
parameters, the classical work of Busse [1^ on the stability of various convective 
states in a infinite fluid predicts that the stable pattern is a set of parallel rolls. 



In the numerical calculations we consider a circular cell of radius R = 327r, which 
corresponds to an aspect ratio F = R/ir = 32. A square grid with A^^ nodes has been 
used with spacing Aa; = Ay = 647r/N, and A^ = 512. We approximate the boundary 
conditions on by taking il){r.,t) = for ||r|| > R, where r is the location of a 
node with respect to the center of the domain of integration. In order to study the 
formation of the chaotic spiral pattern from the conducting state, we choose as initial 
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condition 'il){r,t = 0) a random variable, Gaussianly distributed with zero mean and 
a variance 0.001. 



Our main result is that this model exhibits a spatio-temporal chaotic spiral pattern 
which is remarkably similar to that observed experimentally , (there the system was 
quenched from the conduction state as well). Figure 1 shows a typical configuration 
that we have obtained exhibiting spiral chaos. Dark regions correspond to hot rising 
fluid and white regions to cold descending fluid. Initially, after the control parameter 
is quenched into the regime in which spiral chaos arises, the randomness in the initial 
configuration is rapidly lost. On a time scale t 600, the system self-organizes into 
a regular structure comprising locally rotating spirals that fill the entire cell. 



The two dimensional power spectrum < P{k) >= (\ijj{k)tjj{—k)\'^^ is shown in 
Fig. 2(a), where ip{k) is the spatial Fourier transform of ip{f,t) and <> denotes a 
time average over one horizontal diffusion time |jl9[. The most interesting feature is 



that the intensity of the spectrum appears to be isotropic. Figure 2(b) shows the 
circularly averaged power spectrum P{k). We see that P{k) is broad, skewed and 
peaked at a wavenumber kmax < kc = I. In order to estimate the width of P{k) we 
have fit it to a function of the form A/ (^1 + {k'^ — A^q)^) • The fit obtained is shown 
in Fig. 2(b) as well, and corresponds to A ^ 7.0, ^ ~ 2.4 and ko ^ 0.8. The function 
P{k) is more sensitive to position correlations than orientation correlations in the 
pattern. The dependence of a correlation function that would measure orientation 



correlations on distance need not be the same as that of C(r) [ITS 



We finally show in Fig. 3 the field ( that corresponds to the configuration shown 
in Fig. 1. White and dark regions correspond to clockwise and counterclockwise 
rotation respectively. The Fourier transform of this field shows, as expected, that it 
is peaked around k = and is isotropic. 



Our numerical investigation indicates that both large scale mean flow and large 
aspect ratio play a crucial role in the spontaneous formation of a spiral chaos state. 
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In the absence of the mean flow field, we do not observe it. To study how the 
spiral chaos state depends on the size of the system, a run was conducted in a small 
cell of aspect ratio 16. We observed a globally ordered pattern consisting of one 
two armed spiral, which demonstrates that a large aspect ratio is critical to the 
existence of spiral chaos. We have also studied the role of non-Boussinesq effects 
on the formation of the spiral chaos state. If the non-Boussinesq coupling constant 
g2 = 0, and we start with the same random initial condition and use the same 
parameters (e' = 0.7, Qm = 50, Pr = 1.0 and = 2.0), we observe a similar spiral 
chaos state, demonstrating that the term that models non-Boussinesq effects is not 
necessary for spiral chaos. We have also carried out additional calculations to study 
the effect of the Prandtl number. Starting with exactly the same initial conditions 
and parameters as in Fig.l, but with a larger value of Prandtl number {Pr = 6), 
no spiral chaos pattern is observed, showing that low Prandtl number is essential for 
spiral chaos. The resulting pattern has a labyrinthine aspect. Finally, in order to 
study how generic the spiral chaos state is, we have studied the case in which the 
system was ramped slowly from the conduction state to the convective state. We 
used exactly the same initial conditions and the same parameters as in Fig .1, but 
increased e very slowly up to e = 0.7: e = lO^H for < t < 700, and e = 0.7 for 
t > 700. A similar spiral chaos is observed to form as time increase, although rolls 
tend to persist near the cell boundary. These results suggest that only low Prandtl 
number and large aspect ratio are relevant for the existence of spiral chaos. 
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Figure captions 

Figure 1. Typical configuration of a spiral chaos pattern. The field ip is shown. 
Dark regions correspond to ip > and light regions to ip < 0. The configuration 
shown has evolved from random initial conditions in a cylindrical cell with aspect 
ratio r = 32. The values of the parameters used are g2 — 0.35, — 50 and = 2. 
The parameter e is quenched from to 0.7 in the calculation. The configuration 
shown is at t = 900 pa th- 

Figure 2. (a) Power spectrum in /c-space. The intensity of the spectrum is ap- 
proximately isotropic, (b) Circularly averaged power spectrum. The wavenumber k 
is in units of kc, where kc — 1 is the critical wavenumber. The solid line is the fit 
discussed in the text. 

Figure 3. Vorticity potential field C corresponding to the configuration shown 
in Fig. 1. White and dark regions correspond to clockwise and counterclockwise 
rotation, respectively. 
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